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Abstract 

This paper is divided into two parts. The first part concerns several standard scenarios for 
how short-range spin glasses might behave at low temperature. Earlier theorems of the authors 
are reviewed, and some new results presented, including a proof that, in a thermodynamic sys- 
tem exhibiting infinitely many pure states and with the property (such as in replica-symmetry- 
breaking scenarios) that mixtures of these states manifest themselves in large finite volumes, 
there must be an uncountable infinity of states. 

In the second part of the paper, we offer some conjectures and speculations on possible 
unusual scenarios for the low-temperature phase of finite-range spin glasses in various dimen- 
sions. We include a discussion of the possibility of a phase transition without broken spin-flip 
symmetry, and provide an argument suggesting that in low dimensions such a possibility may 
occur. The argument is based on a new proof of Fortuin-Kasteleyn random cluster percolation 
at nonzero temperatures in dimensions as low as two. A second speculation considers the pos- 
sibility, in analogy to certain phenomena in Anderson localization theory, of a much stronger 
type of chaotic temperature dependence than has previously been discussed: one in which the 
actual state space structure, and not just the correlations, vary chaotically with temperature. 
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1 Introduction 



In this paper we consider Ising spin glass models, in particular the infinite-range Sherrington- 
Kirkpatrick (SK) model [ 1 1 and the nearest neighbor Edwards-Anderson (EA) model on Z d . 
The SK Hamiltonian (for N spins) is 

H N = -(1/VN) E (!) 

\<i<j<N 

where the couplings Jy are independent, identically distributed random variables chosen, e.g., 
from a Gaussian distribution with zero mean and variance one. 

In a series of papers, Parisi and collaborators @J |6j] proposed, and worked out the con- 
sequences of, an extraordinary ansatz for the nature of this phase. Following the mathematical 
procedures underlying the solution, it came to be known as replica symmetry breaking (RSB). The 
starting point of the Parisi solution was that the low-temperature spin glass phase comprised not just 
a single spin-reversed pair of states, but rather "infinitely many pure thermodynamic states" H, 
not related by any simple symmetry transformations. 

What is primarily of interest is the distribution of the overlaps between two SK 'pure states' 
a and (3, defined as 

1 N 

9aP = Z)<0i)a(Oi)/9 , (2) 

iV 8=1 

where (-) a is a thermal average in pure state a, and dependence on the coupling realization J and 
temperature T has been suppressed. The overlap distribution is constructed by choosing, at fixed 
N and T, two of the pure states a and (3 present with probabilites W a and Wp in the Gibbs state. 
The probability that their overlap lies between q and q + dq is then given by the quantity Pj(q)dq, 
where 

Pj(Q) = E E W$W$5(q - q a p) • (3) 

a /3 

As before, we suppress the dependence on T and N for ease of notation. The average P(q) of 
Pj(q) over the disorder distribution is commonly referred to as the Parisi overlap distribution, and 
serves as an order parameter for the SK model. 

The EA Hamiltonian in zero external field is given by 

H — — E JxyOxOy , (4) 

\x— y\=l 

where the nearest-neighbor couplings J xy are defined in exactly the same way as the in the 
SK Hamiltonian dU). In this model, and unlike in the SK model, thermodynamic pure, mixed, 
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and ground states are standard, well-defined (see, e.g., JEHU) objects, constructed according to 
well-established prescriptions of statistical mechanics lfT0lfTTlfT^fT3lfT4lfT5lfT6ll . 

Now suppose — as has often been conjectured for finite-range spin glasses ifTTl ll"8ll — that 
there are many pure states not simply related to each other by symmetry transformations. If this is 
the case, it was shown in [19| that local variables, such as correlations, will vary chaotically and 
unpredictably as volume size changes (with, say, periodic boundary conditions on each volume). 
Is it then even possible to describe the nature of large finite volume systems via a single infinite 
volume object, and if so, how? It turns out that an object, called the metastate [20| by the authors, 
can be constructed to capture the nature of the behavior inside finite systems as volume increases. 
The metastate is a useful tool that describes the empirical distribution of local variables as volume 
increases without bound, and is described in more detail in the companion paper lETll and in ll8ll2"Ul 
1221 123 I2H 123. The seemingly chaotic behavior with increasing volume is modelled by random 
sampling from the metastate, regarded as a probability measure on the space of Gibbs states. Using 
the metastate approach, we now examine more closely the nature of the low-temperature spin glass 
phase. 

2 The trinity of scenarios 

There are three scenarios that have received the major share of attention in the current literature. 
Here we refer only to those dealing with the pure state structure of the low-temperature spin glass 
phase. Other pictures have also received a great deal of attention, particularly the excited state 
scenario of Krzakala-Martin [26| and Palassini- Young [27 1. Proving or disproving these or related 
pictures is important for achieving a thorough understanding of the low-temperature physics of 
spin glasses. However, because it can be shown [28, 29 1 that these pictures describe excitations 
that do not alter pure state structure, their presence or absence can be logically incorporated into 
any of the three pictures that we are about to describe (although on heuristic grounds they appear 
incompatible with the RSB picture). 

Low-temperature pictures of spin glass long-range order and broken symmetry start with an 
assumption about the number of pure states Af({3, d) (which is the same for almost every coupling 
realization J — see, e.g., lE2l ~). They assume also a putative critical (inverse) temperature f3 c (d) < 
oo separating a paramagnetic phase for (3 < f3 c from a spin glass phase for (3 > f3 c . Although 
there is good numerical IfTTl |3*U1 |3"T1 l3*2l and some analytical [[331 13~31 evidence that above some 
lower critical dimension d\ there does exist such a finite f3 c (d), there is as yet no proof or even a 
strong physical argument supporting such a conjecture. Moreover, there is no logical reason why 
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there cannot be two or more phase transitions in some dimensions. However, we will not attempt 
to enumerate all possibilities here; we will confine ourselves to the most likely scenarios. We 
defer to SecUJa consideration (mostly for the fun of it) of some of the more outlandish sounding 
possibilities. 

The actual value of Af(fi, d) is not rigorously known at large (3 for any d > 2. There does 
exist a rigorous argument ll35l supporting — but not completely proving — the conjecture that 
jV(oo, 2) = 2, and a heuristic argument tl36l supporting the conjecture that J\f(f3 < oo, 2 < 
d < 8) < 2. It is generally assumed (but also not proved) that spin-flip symmetry is broken for 
f3 c < (3 < oo, so that pure states come in global spin-flip reversed pairs. 

Given all these assumptions, pure state scenarios generally assume either a single pair of pure 
states or infinitely many. Again, there is at this time no argument proving that one cannot have, 
say, JV(P, d) = 20; it's just difficult to imagine why such a scenario should occur. There is, 
however, a reasonably strong heuristic argument [24 1 indicating that if JV(P, d) = oo, it must be 
an uncountable infinity (we will prove this in Sec.|3]for the RSB picture). 

While an assumption about JV(P, d) is the starting point for each of the three pictures we now 
describe, they are each much more than a simple assertion about the number of pure states. In 
particular, there are potentially many 2-state or many-state pictures that correspond to none of 
these three scenarios. Here the only aspect of the three on which we will focus is the relationship 
among the pure states. 

2.1 Heuristic Description 

The scenarios are, in order of increasing complexity, the scaling/droplet (SD) picture ll37l l38l l39l 
gOlSD |42l|43l, the chaotic pairs (CP) picture 00[I2EQl|23l24l, and the RSB picture EE 00 
[TBI . The first is a 2-state scenario, while the second and third are many-state scenarios. The pure 
state structure in these pictures is normally described through the Parisi overlap function P(q) lfl8l . 
This function needs to be used with great care in describing pure state structure in short-range 
models; see ED for a discussion of some of the pitfalls and problems that can occur in its 

applications. 

In the next subsection we provide precise definitions of the pure state structure of these pictures; 
we limit ourselves here to brief heuristic descriptions. As already noted, scaling/droplet is a two- 
state picture. The overlap distribution is therefore simply a pair of ^-functions at q = ±qEA, where 
qEA is the Edwards- Anderson order parameter [2J, regardless of whether replicas and overlaps 
are taken before or after the thermodynamic limit is taken, as shown in Fig. 1 (see below and 
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Figure 1: (From Ref. J9|.) The spin overlap function P(q) at T < T c for: (a) a two-state picture 
when replicas are taken before the thermodynamic limit; (b) the many-state chaotic pairs picture 
when replicas are taken before the thermodynamic limit; (c) a two-state picture when replicas are 
taken after the thermodynamic limit; (d) the many-state chaotic pairs picture when replicas are 
taken after the thermodynamic limit (conjectured). 

especially j8l|9l[22J for a more detailed discussion of the difference between the two procedures). 

The CP picture considers the possibility of an infinite number of pure state pairs, but with a 
trivial overlap structure: all large, finite volumes would display an overlap structure equivalent to 
the SD picture because only a single pair of pure states appears in each of these volumes. The 
actual pure state pair, however, varies chaotically as volume changes because different volumes 
select different members from the infinite ensemble of such pairs. The overlap distribution for the 
infinite volume in the CP picture, constructed by choosing replicas and taking their overlaps after 
going to the thermodynamic limit, would then presumably be a simple 5-function at zero overlap 
(cf. Fig. 1). 

There are actually two RSB pictures, which we have called 'standard' and 'nonstandard' 
SK ll20l l23l l24l l44l . The first chooses replicas after taking the thermodynamic limit, and the 
second before. Both have nontrivial overlap structures, which will be described in the next subsec- 
tion. 

Before turning to that, we need to say a few words about the dimension dependence of the three 
pictures. The only one of these with specific predictions is the SD picture, which asserts that in 
every finite dimension where spin-flip symmetry is broken, there is only a single pure state pair. 
The CP picture does not make any corresponding claim; it merely asserts that if Af(/3, d) = 00 at 
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some (/3, d), then the overlap structure must be trivial in the manner specified above. If the RSB 
picture correctly describes the EA model at d — oo, then both the SD and CP models agree that 
the upper critical dimension = oo for the EA model. 

The RSB picture is slightly less vague about its dimension-dependence; it does apparently 
assume that there exists a strictly finite d u c for the EA model, and that replica symmetry is broken 
in the nontrivial manner it specifies for all d > d". The precise value of d u c within this picture 
remains uncertain, but it seems to be higher than two and less than or equal to three P31B5I . 

2.2 Description via the Metastate 

We now turn to a precise description of the three competing pictures, which is greatly facilitated 
by using the metastate. 

In SD, there is only a single pair of pure states, and these are the same in every large volume; 
the overlap distribution function Pj(q) in a volume A L therefore simply approximates a sum of 
two 5-functions at ±qEA, as shown in Fig. 1(a). In the CP picture, each finite-volume Gibbs state 
Pj will still be approximately a mixture of a single pair of spin-flip-related pure states, but now 
the pure state pair will vary chaotically with L. Then for each A L , Pj(q) will again approximate 
a sum of two ^-functions at ±qEA- 

So chaotic pairs resembles the scaling/droplet picture in finite volumes, but has a very different 
thermodynamic structure. It is a many-state picture, but differs from the RSB picture (see below) 
in that only a single pair of spin-reversed pure states pj L , p °j , appears in a large volume with 
symmetric boundary conditions, such as periodic. In other words, for large L, one finds that 

Pf ™\pj L + \$F > (5) 

and the pure state pair (of the infinitely many present) appearing in a particular finite volume A t 
depends chaotically on L. That is, the periodic b.c. metastate is dispersed over many T's (in fact, 
an uncountable infinity, as we shall see below; this also occurs in the RSB picture), but (unlike in 
RSB) each T is a trivial mixture of the form Y = T a = \p c j + \p ( j- The overlap distribution for 
each r is the same: P r = \8{q — q^A) + \&{q + Qea)- It is interesting to note that there is a spin 
glass model (the 'highly disordered model' J47l|48l|49|) that appears to display just this behavior 
in its ground state structure above eight dimensions. 

We now turn to the standard and nonstandard mean-field-like scenarios. The standard SK pic- 
ture is perhaps most concisely described in the introduction of [ 50]. It requires a Gibbs equilibrium 
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measure pj{p) which is decomposable into many pure states Pj(cr): 

Pj(a)=J2W^p a j(a). (6) 

a 

In this picture replicas are taken after going to the thermodynamic limit. That is, one chooses 
cr and a' from the product distribution pj(cr)pj(a'), and then the overlap can be defined as 

Q= lim \Al\- 1 ]T a x a' xJ (7) 
l— >oo — r 

x€A L 

where |A^| is the volume of the cube A^. 

Suppose that a is drawn from pure state pj and a' from pj. Then © equals its thermal 
mean P4l 

gj 7 = ^lmi^ |A L | _1 ^ (o- 1 .) a (a z ) 7 , (8) 
and so the overlap distribution Pj(q) is given by 

a, 7 

According to this picture, the Wj's and q'j's are non-self-averaging quantities, except when 
a = 7 or its global flip, where gj 7 = ±qEA- The average -P(g) of Pj{q) over the disorder 
distribution v of the couplings is then a mixture of two delta-function components at ±qEA and a 
continuous part between them. 

There is a technical problem (caused by chaotic size dependence |[T9lO in the construction of 
Pj(cr) that can be overcome by using the periodic b.c. metastate ft™ (in fact, any coupling- 
independent metastate would do). One can construct P4ll a state pj(cr) which is the average over 
the metastate: 

Pj (a) = ( T{a)Kj{Y) dT . (10) 



One can also think of this pj as the average thermodynamic state, N~ 1 (p^ + p { j 2) + ..., p { j N) ), 
in the limit — > oo. It can be proved G21E51 that pj(a) is a Gibbs state. 

Numerically one constructs overlaps without constructing Gibbs states at all. Such a construc- 
tion (similar to that above) is described in [ 44l , and leads ultimately to the same conclusion. 

But this picture can be rigorously ruled out for the EA model [44|, as will be shown in Sec.|3] 
So the most natural (and usual) interpretation of a mean-field-like picture cannot be applied to 
short-range spin glasses. The question then becomes: are there alternative, less straightforward 
interpretations? To address this question, we constructed in EDI (see also |9j |22l |23l I24l0 
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q EA q q EA q EA q 

Figure 2: (From Ref. [8].) The overlap distribution, at fixed J , in two different volumes Ai and 
A 2 in the nonstandard SK picture. 

an alternative mean-field-like picture, the 'nonstandard SK model', which we hereafter refer to 
simply as the RSB picture. This picture clarifies as well how broken replica symmetry should be 
interpreted for the SK model. 

Consider again the PBC metastate, although, as always, almost any other coupling-independent 
metastate will suffice. The RSB picture assumes that in each volume A^, the finite- volume Gibbs 
state pj t Li is well approximated deep in the interior by a mixed thermodynamic state Y^ Li \ de- 
composable into many pure states p aL : 

r w = E^ (id 

where explicit dependence on J is suppressed. 

As in the chaotic pairs picture, the mixed states change in some "chaotic" fashion with 
Li. Furthermore, each mixed state is presumed to have a nontrivial overlap distribution 

PtW= E K%Ki^(Q ~ Qa H p Li ) (12) 

of the form shown in Fig. 2. Moreover, the distances among any three pure states within a partic- 
ular T are assumed to be ultrametric lffl~8l . 

We conclude with a brief note about zero temperature. In each A L (with periodic b.c.'s) there 
can only be a single ground state pair at T = (because we are considering Gaussian couplings 
rather than a ± J model). If scaling/droplet holds, then this pair is the same for all large L; if 
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infinitely many ground state pairs exist, then the pair changes chaotically with L. In this respect 
the behavior of CP is the same at both zero and nonzero temperatures (below T c ), and the same 
can be said for SD. So the overlap functions of CP and RSB differ only at positive temperature: 
the mean-field RSB picture at T > has the r^ L ^ in each volume exhibiting a nontrivial mixture of 
pure state pairs as in dTTb . while in chaotic pairs the r^ L ^ appearing in any A L consists of a single 
pure state pair, as in ©. There remains, however, an important difference between the CP and 
RSB pictures at zero temperature, in that the space-filling interfaces between ground states should 
have energies that scale differently; see and the companion paper in this volume ET1 for a more 
detailed discussion. 

3 Replica Symmetry Breaking for Short-Range Models 

We begin by eliminating the standard SK picture from further consideration; the detailed rigor- 
ous proof appears in ll44ll . Using the torus-translation symmetry of the periodic b.c.'s, one can 
show that the Gibbs state pj(cr) is tramlation-covariant; that is, pja(a) = pj(a~ a ), or in terms 
of correlations, (a x ) ja = (a x ^ a ) j, and similarly for n-point correlations. Translation covari- 
ance of pj immediately implies, via Eqs. Q-©, translation invariance of Pj. But, given the 
translation-ergodicity of the underlying disorder distribution u, it immediately follows that Pj(q) 
is self-averaging, and equals its distribution average P(q) for a.e. J . The impossibility of a 
non-self-averaging Pj(q) can also be shown for other coupling-independent b.c.'s, where torus- 
translation symmetry is absent, using methods described in ll35ll . 

This leaves nonstandard SK as a maximally allowed mean-field-type picture. Before discussing 
its viability, we explain why neither CP nor nonstandard SK remains viable if modified to have 
only a countable infinity of pure states. We do this first by a heuristic argument valid for both 
CP and nonstandard SK, and then by a rigorous argument valid only for nonstandard SK. The 
conclusion of this exercise is that both alternatives to the SD scenario of a single pure state pair 
require uncountably many (i.e., a continuum of) pure state pairs for a single fixed J . 

Consider p j, the average over the periodic b.c. metastate k j. Changing from periodic boundary 
conditions to antiperiodic, or to any of the many other 'partially antiperiodic' boundary conditions 
related to periodic ones by a gauge transformation, leaves k j and consequently p j unchanged ll24ll . 
On heuristic grounds, unless pj is of the SD form with only a single pair of pure states, then this 
lack of dependence on boundary conditions should hold only if pj is a uniform mixture over 
infinitely many pure states pj — i.e., the relative weights of all pj's in pj are equal. But that is 
clearly impossible if the number of p°j% is countably infinite. 
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In the case of an SK-type picture, where the T's appearing in the metastate k j involve nontrivial 
weights Wj, one can go further and prove that nontrivial Wj's require uncountably many p^'s to 
appear in kj. The argument proceeds as follows. If there are only countably many pj's, then 
the corresponding weights (say, ordered by magnitude) in pj will be measurable and translation- 
invariant functions of J . By translation-ergodicity, this means that they are in fact independent of 
J and so remain the same if finitely many couplings J xy change by amounts AJ xy . Moreover, one 
can consider the finitely many a's corresponding to, say, the largest k weights in pj, and then the 
distribution of their weights within the T's of the metastate kj is also independent of J . 

On the other hand, this leads to a contradiction, because a Y in the nonstandard SK metas- 
tate (see Eq. (fTTT) ') is a nontrivial mixture of pure states and their weights. By the Aizenman- 
Wehr transformation ll2"51l . a pure state p a transforms to a pure state p a > under such a finite change 
J — > J + A J, and its weight W a within T correspondingly changes according to 

Wa - W a , = r a Wj r,W, (13) 

7 

where 

r a = (exp(/3 AJxyVxVy))a ■ (14) 

(xy) 

For pure states a and 5 that are not identical one can find a choice of A<J, with all AJ xy 
small, such that r a ^ r$ . Then for all T's in which W^r) and W^r) are nonzero, the ratio 
W a >(T)/Ws'(T) is changed from W a (T)/Ws(T) by the same factor r a /r s , that is close to I. But 
this implies that J — >• J + A J will change the distribution of the weights within the metastate, 
leading to a contradiction. We have thus sketched the proof of: 

Theorem 1. The PBC metastate k j for the EA spin glass cannot assign strictly positive prob- 
ability to T's whose pure state decompositions satisfy both: 

(i) r = E« W%p a j, with not all Wj = 1/2, and 

(ii) over all these T's only countably many pj's appear. 

An immediate consequence of Theorem 1, when combined with our earlier arguments, is that 
the only remaining mean-field-like picture is the nonstandard SK model with a (periodic b.c.) 
metastate k j whose average pj is supported on an uncountable infinity of pure states. 

We cannot at this point rule this scenario out rigorously, but can provide a strong heuristic 
argument, based on a rigorous result, that makes it very unlikely. The rigorous result is the already 
mentioned invariance of the metastate (24|: two metastates constructed using gauge-related b.c.'s 
(e.g., periodic and antiperiodic, or any two randomly chosen fixed b.c.'s) are identical. This makes 
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it difficult to see how any many-state picture can have a pj supported on anything other than a 
uniform distribution of pure states. But if this did occur for a particular J , an Aizenman-Wehr 
transformation suggests that the uniformity would be destroyed for a finitely different J' . (For 
detailed arguments, see f2M .) 

Why doesn't this argument also rule out chaotic pairs? Because in the chaotic pairs picture, as 
in scaling/droplet, there are in each r( L ) only two pure states (although in CP the pair depends on 
L), each with weight 1/2. All even correlations are the same in any pair of flip-related pure states, 
so, by Eqs. (fT3t and (fl4ll . any change in couplings leaves the weights unchanged. 

We conclude that the only viable many-state scenario for short-range Ising spin glasses is the 
CP picture with an uncountable infinity of pure states. This picture has trivial replica symmetry 
breaking (cf. Fig. 1) and consequently a very simple overlap structure. 

4 Wild Possibilities 

In Sec. El we presented the most likely scenarios for the pure state structure of the spin glass phase 
in short-range models. One of these, the RSB picture, although a priori plausible, especially given 
its relevance to the SK model, was rigorously ruled out in two of its possible versions and excluded 
heuristically in its final remaining version. However, there are other, more exotic scenarios that 
might also occur but do not seem to have been considered. Here, for the sake of completeness 
(and perhaps a little whimsy) we suggest two. In order to keep the list of possibilities relatively 
constrained, we will assume in all that: a) a thermodynamic phase transition does exist above 
some lower critical dimension d l c , and b) in a fixed dimension the low-temperature phase does not 
alternate among different scenarios, such as SD and CP, as temperature is lowered (although the 
third and most exotic scenario may be regarded as a kind of phase of extreme alternation). 

4.1 Phase transition without broken spin-flip symmetry 

As mentioned in Sec.|2j despite decades of effort and considerable numerical support lfT7ll30ll3Tl 
l32l . there remains no proof of a phase transition in EA spin glasses. In this section we provide a 
proof that shows that a necessary (though not sufficient) condition for broken spin-flip symmetry 
at T > is satisfied in lattices with dimension as low as two (e.g., a triangular lattice). We 
also discuss a corresponding sufficient condition. Of course, this (minimal) symmetry-breaking is 
assumed in all scenarios discussed in Sec. but here we speculate on the possibility of a phase 
transition without the appearance of multiple Gibbs states, and discuss why the presence of the 
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necessary condition (for broken spin-flip symmetry) at some (/?, d) where the sufficient condition 
is absent suggests just such a possibility. 

Although several approaches are possible, we use here the Fortuin-Kasteleyn random cluster 
(RC) representation EJEI extended to non-ferromagnetic models (see, e.g., 11531 ). which relates 
the statistical mechanics of Ising (or Potts) models to a dependent percolation problem. We let E d 
denote the set of bonds (x, y) in some <i-dimensional lattice L d , each with corresponding coupling 
J xy . The RC approach introduces parameters p xy G [0, 1) by: 

p xy = 1- exp[-/3| J xy \) . (15) 

The RC distribution is a probability measure /xrc on {0, l} Ed , that is, on 0- or 1 -valued bond 
occupation variables n xy . It is one of two marginal distributions (the other being the ordinary 
Gibbs distribution) of a joint distribution on Q = { — 1, +l} Ld x {0, l} Ed of the spins and bonds 
together, and is given (formally, in the infinite system) by 

//rcKv) = Z R h 2 # «"-» fx ind ({n xy }) lu{{n xy }) , (16) 

where Z RC is a normalization constant, 4^{{ n xy]) is the number of clusters determined by the 
realization {n xy }, f^md({ n xy}) is the Bernoulli product measure corresponding to independent oc- 
cupation variables with ^i n d({ n xy — 1}) = Pxy, an d It/ is the indicator function on the event U in 
{0, l} E<i that there exists a choice of the spins a x so that J xy n xy a x a y > for all (x, y) Il54ll55ll56ll . 
U is the event that there is no frustration in the occupied bond configuration. Finite-volume ver- 
sions of the above formulas, with specified boundary conditions, are similarly constructed. 

The mapping of this formalism to ferromagnets follows from formulae (that do not hold for 
nonferromagnets) such as ll57l 

(vxcr y ) = Vrc(x y) , (17) 

where (cr x a y ) is the usual Gibbs two-point correlation function and /iRc(£ <-> y) is the RC prob- 
ability that x and y are in the same cluster. Similarly, with appropriate RC (wired) boundary 
conditions used for fi RC , one has 

(o x ) + = fi RC (x <-> oo) . (18) 

It follows that, for ferromagnets, a phase transition from a unique (paramagnetic) phase at low 
P to multiple infinite volume Gibbs distributions at large (3 is equivalent to a percolation phase 
transition for the corresponding RC measure. 
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For spin glasses (or other nonferromagnets) the situation is less straightforward. Now for two 
sites x and y, (flTt becomes 

(<7x<7y) = (lx~aoV(B,y))RC] Vi^^V) = II Sgn(J^ v ') , (19) 

(x'y')eC 

where C is any path of occupied bonds from x to y. By the definition of U, any two paths C and C 

in the same cluster will satisfy Y[(x'y>)&c s S n (Jx'y') = H(x'y')eC' s S n (Jx' y ')- 

It is no longer evident that RC percolation is sufficient to prove broken spin-flip symmetry. 
Consider the case of a finite volume A L with fixed boundary conditions, i.e., a specification a x = 
±1 for each a x G 8Al- For the ferromagnet, by first choosing all a x = +1 and then all W x = — 1, 
one can change the sign of the spin cr at the origin even as L — > oo. That is, boundary conditions 
infinitely far away affect a , which is a signature of the existence of multiple Gibbs distributions. 

But for the EA spin glass, it is not clear whether, even in the presence of RC percolation, there 
exist any two sets of boundary conditions with the same effect. Although the infinite cluster in any 
one RC realization is unique, different RC realizations can have different paths from «-> 8Al, 
leading to different signs for a . So percolation might still allow for (cr ) — ► as L — ► oo. 

However, it is easy to see that single RC percolation is at least a necessary condition for mul- 
tiple (symmetry-broken) Gibbs phases. The contribution to the expectation of a from any finite 
RC cluster is zero: if a spin configuration a is consistent with a given RC bond realization within 
such a cluster, —o is also consistent and equally likely. As a consequence, at a (f3, d) where RC 
percolation does not hold, (a ) = in infinite volume. We note that a slightly stronger version 
of this argument [56| proves that the transition temperature for an EA ±J (or other) spin glass, 
if it exists, is bounded from above by the transition temperature in the corresponding (disordered) 
ferromagnet. 

Is there a modification of this approach that could lead to a proof of multiple Gibbs phases? One 
such possibility is what might be called double RC percolation. Here one expands the sample space 
Vt to include two independent copies of the bond occupation variables (for a given J configuration), 
and defines the variable r xy = n xy n' xy , where n xy and n' xy are taken from the two copies. One 
then replaces percolation of {n xy } in the single RC case with percolation of {r xy }. It is not hard 
to see that this would be a sufficient condition for the existence of multiple Gibbs phases (and 
consequently, for a phase transition). 

Single RC percolation for the EA ± J Ising model in d > 2 11551 has been proved l!53l . Here we 
sketch the outline of a proof (simpler than the one in (521) showing that one has single RC perco- 
lation even in d = 2 — e.g., on the triangular lattice. This is already interesting for the following 
reason: because we are not aware of strong numerical evidence of multiple Gibbs states for this 
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geometry, this may be considered evidence that single RC percolation can occur in spin glasses 
without broken spin-flip symmetry. But there is an even more interesting potential consequence, 
which we will discuss below. 

The proof uses a standard Fortuin-Kasteleyn-Ginibre (FKG) argument ll58ll using correlation 
inequalities. Let E L be the edge set confined entirely within a volume A^. It can be shown that 
there exist probability measures v L on {0, l} Ei such that, if f({n}) and g({n}) are nondecreasing 
real functions (that is, they do not decrease when any {n} — > {n'}, where every n xy = 1 corre- 
sponds to n' xy = 1, but n xy = can correspond to either n' = or 1), then they are positively 
correlated: 

{fg) VL >{f)uM^. (20) 

One measure satisfying this property is the independent product measure [ 59 ] . Consequently, 
the marginal distribution (i.e., averaged over the coupling (signs) in the ± J model) of satisfied 
bonds (i.e., using the "satisfaction" variables n xy = 1 (or else 0) if J xy a x a y > 0) at (3 = 0, 
satisfies this property. Operationally, one can think of constructing this set by choosing the spins 
first through independent flips of a fair coin, and then choosing the sign of each bond in the same 
way. One then has independent, density- 1/2 bond occupation, which percolates on the triangular 
lattice EDI ED- 

Consider now the satisfaction variables n b for (3 > 0. It is not hard to show that, at any such 
fixed [3, and for any bond b = (xy), 

P(n b = 1| {ny : b' ^ b}) > 1/2 + e(/3) , (21) 

where e(/3) is a nonnegative function of [3. Eq. (OTT) implies that percolation of satisfied bonds 
at finite temperature dominates (in the FKG sense) percolation at infinite temperature. Since the 
RC variables n b are obtained from the satisfaction variables n b by a slight (for large (3) dilution, it 
follows that at sufficiently low (but nonzero) temperature one has single RC percolation. 

It is presumably the case that on the triangular lattice there is no broken spin-flip symmetry, 
and also only a single Gibbs state, at all nonzero temperatures. But it is worth entertaining the 
possibility that single, but no double RC percolation, does imply some sort of phase transition but 
with a single Gibbs state at all nonzero temperatures. 

Let [3 C < oo be the inverse RC percolation transition temperature for the EA model on the 
triangular lattice. Consider again the expectation of the spin at the origin, (cr ), in a volume A L 
with plus boundary conditions. For (3 < (3 C , the probability of the site belonging to a cluster 
reaching the boundary is bounded from above by co{(3)e~ Cl ^ L , where each Ci{(3) > is a finite 
constant. Therefore (a ) < co((3)e~ cl ^ L . 
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For (3 > (3 C , however, the probability that the origin belongs to a cluster reaching the boundary 
is bounded away from zero as L — > oo. In order for no phase transition to take place, there 
must be at fixed [3 an almost exact cancellation between those RC realizations of 'positive' (in 
the product of the signs of the couplings) and 'negative' paths from the origin to the boundary. It 
is at least conceivable that, while a — > as L — > oo, it falls off slower than exponentially — 
perhaps as a power law in L. This would imply a phase transition from a paramagnetic phase at 
high temperature to a phase at low temperature with a unique Gibbs state, but one where two-point 
correlations decay as a power law: (cr x a y ) ~ \x — y\~ a W as \x—y\ — > oo. This would be analogous 
to the Kosterlitz-Thouless phase transition for 2D XY models ll62l . 

4.2 Highly chaotic temperature dependence 

We discussed earlier [19| the correspondence between multiple Gibbs state pairs and the appear- 
ance of chaotic size dependence of correlations in an infinite sequence of volumes with, say, pe- 
riodic boundary conditions. A possibly related phenomenon is the speculated presence in spin 
glasses (both infinite- and short-range) of chaotic temperature dependence (CTD). Roughly speak- 
ing, this refers to an erratic behavior of correlations, upon changing temperature, on lengthscales 
that diverge as the temperature increment goes to zero. CTD was predicted [39||43l|63l for the 
EA spin glass within the SD context; it also seems to be implied by the RSB theory Il64ll65ll66l. 
More recent numerical and analytical work (see ll67l and references therein) have led to claims 
that CTD is not present in either the SK or the EA model (see also [68|), although [69| allows 
for the possibility of a weak effect at large lattice sizes. A perturbative approach 11701 observes a 
very small effect at ninth order. At this time the issue remains unresolved. Its potential presence in 
spin glasses is interesting, however, and represents a qualitatively new thermodynamic feature of 
at least some types of disordered systems. 

In this section we raise the possibility of a far stronger version of CTD, which we will call 
highly chaotic temperature dependence (HCTD) [36|. Unlike 'ordinary' CTD, where correlations 
behave in a chaotic fashion as temperature changes but the global pure state structure does not, in 
HCTD the number of pure states in the periodic b.c. metastate does not behave in a continuous or 
even monotonic fashion as temperature is lowered. This picture departs radically from any other 
that has appeared so far in the literature. 

The HCTD scenario is summarized as follows. As in all other pictures, there exists a determin- 
istic T c for a.e. J, such that for all T > T c there is a unique (paramagnetic) Gibbs state. Suppose 
one were now to choose an arbitrary (nonzero) T < T c . Then with probability one (i.e., for almost 
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every J), there is again a unique infinite volume Gibbs state (though maybe not with exponential 
decay of truncated correlations). Nevertheless, as temperature is lowered from T c to for a fixed 
typical J, there would be a (countably infinite) dense set of temperatures, depending on J, with 
broken symmetric pure phases for that J at those temperatures. 

Now this scenario might be dismissed, with some justification, as the authors' fevered imag- 
inings, but it should be noted that just this sort of phenomenon does happen in other disordered 
systems. Consider Anderson localization in one dimension, with the random Schrodinger opera- 
tor CD 

H w = -V 2 + V (x) + XV u (x) , (22) 
where A > 0, V (x) is a bounded potential periodic in x, and 

oo 

V u {x)= £ r%U(x-x n ). (23) 

n=— oo 

In (l23l n runs over the integers, U(x — x n ) is a nonegative localized potential centered at lattice 
site x n , and the r\ n are i.i.d. random variables uniformly distributed in [0, 1]. 

Operators of the type given by d22l are known to have certain properties - see, e.g. ll72ll . If one 
first chooses a specific energy, then for a.e. realization V w of the random potential, that energy is 
part of its continuous spectrum (i.e., is not an eigenvalue of Hu,). On the other hand, if one first 
chooses a specific realization V w , then there is almost surely a (countable) dense set of eigenvalues 
(in some appropriate energy interval) of H^, and of course this set depends on u. 

Returning to the EA model in (say) three dimensions, it is amusing to speculate that a phe- 
nomenon like that described above for localization might occur for EA spin glasses, as follows: a) 
For every arbitrarily chosen T, there is a unique Gibbs state for a.e. J; but b) there exists a (deter- 
ministic) T c , such that for a.e. J , the set of temperatures T such that there are multiple Gibbs states 
(e.g., with q E A ^ 0) for that J is dense in the entire temperature interval [0, T c ). By an application 
of Fubini's theorem (73 1, property a) would necessarily imply that the set of such T(j7")'s would 
have zero Lebesgue measure in the temperature line. 
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